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We discuss the problem of the formation of a large-scale magnetic field in the accretion disks 
around black holes, taking into account the nonuniform vertical structure of the disk. The high 
electrical conductivity of the outer layers of the disk prevents the outward diffusion of the mag- 
netic field. This implies a stationary state with a strong magnetic field in the inner parts of the 
accretion disk close to the black hole, and zero radial velocity at the surface of the disk. Magnetic 
jet collimation is considered, when the jet radius is hold due to magneto-torsional oscillations. 
The range of parameters is found where jet radius is oscillating, regularly or chaotically, within 
restricted values. 
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1. Introduction 

Quasars and AGN contain supermassive black holes, about 10 HMXR contain stellar mass 
black holes - microquasars. Jets are observed in objects with black holes: collimated ejection from 
accretion disks. 

Early work on disk accretion to a black hole argued that a large-scale magnetic field of, for 
example, the interstellar medium would be dragged inward and greatly compressed by the accreting 
plasma [Q, |6[ [J. Subsequently, analytic models of the field advection and diffusion in a turbulent 
disk suggested, that the large-scale field diffuses outward rapidly [(l2|, [RJ, and prevents a significant 
amplification of the external poloidal field by electrical current in the accretion disk. This has 
led to the suggestion that special conditions (non-axisymmetry) are required for the field to be 



advected inward [|14p. The question of the advection/diffusion of a large-scale magnetic field in a 
turbulent plasma accretion disk was reconsidered in [|J], taking into account its nonuniform vertical 
structure. The high electrical conductivity of the surface layers of the disk, where the turbulence is 
suppressed by the radiation flux and the relatively high magnetic field, prevents outward diffusion 
of the magnetic field. This leads to a strong magnetic field in the inner parts of accretion disks 
around black holes. 



2. The fully turbulent model 

There are two limiting accretion disk models which have analytic solutions for a large-scale 
magnetic field structure. The first was constructed in [||, ^| for a stationary non-rotating accretion 
disk. A stationary state in this disk (with a constant mass flux onto a black hole) is maintained by 
the balance between magnetic and gravitational forces, and thermal balance (local) is maintained 
by Ohmic heating and radiative conductivity for an optically thick conditions. The mass flux to 
the black hole in the accretion disk is determined by the finite conductivity of the disk matter and 
the diffusion of matter across the large-scale magnetic field as sketched in Fig.l. The value of the 
large-scale magnetic field in stationary conditions is determined by the accretion disk mass, which 
in turn is determined by the magnetic diffusivity of the matter. It is widely accepted that the laminar 
disk is unstable to different hydrodynamic, magnetohydrodynamic, and plasma instabilities which 
implies that the disk is turbulent. In X-ray binary systems the assumption about turbulent accretion 



disk is necessary for construction of a realistic models [|13j]. The turbulent accretion disks had been 
constructed also for non-rotating models with a large-scale magnetic field. A formula for turbulent 
magnetic diffusivity was derived in [0]. similar to the scaling of the shear a-viscosity in turbulent 



accretion disk in binaries [|13|], where the viscous stress tensor component t r § = aP, with a < 1 a 
dimensionless constant and P the pressure in the disk midplane. disk. Using this representation, 
the expression for the turbulent electrical conductivity a t is written as 

c 2 

Ot = 7=- (2.1) 

aAnh^/P/p 

Here, a = a\a^. The characteristic turbulence scale is £ = a\h, where h is the half-thickness 
of the disk, the characteristic turbulent velocity is v t = OLi^/Pj 'p. The large-scale magnetic field 
threading a turbulent Keplerian disk arises from external electrical currents and currents in the 
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Figure 1: Sketch of the poloidal magnetic field threading an accretion disk. The field strength increases 
with decreasing radius owing to flux freezing in the accreting disk matter, from [0] . 



accretion disk. The field generated by the currents in the disk can be much larger than that due 
to the external currents. The magnetic field may become dynamically important, influencing the 
accretion disk structure and leading to powerful jet formation, if it is strongly amplified during the 
radial inflow of the disk matter. It is possible only when the radial accretion speed of matter in the 
disk is larger than the outward diffusion speed of the poloidal magnetic field due to the turbulent 
diffusivity rj t = c 2 /(4no t ). Estimates in [ jl2[ ] have shown that for a turbulent conductivity 



the outward diffusion speed is larger than the accretion speed. Thus it appears that there is no 
large-scale magnetic field amplification during Keplerian disk accretion. Numerical calculations in 
[ |12| ] are reproduced analytically for the standard accretion disk structure which can be written as 
(e.g.[01) 



M Anpv r rh , h — — , v s 



- , AnrhaP - 
P 



■M(j-ji, 



2 r 



2aT 4 c 
3Kph 



(2.2) 



Far from the inner disk boundary the specific angular momentum is j » jt n . The characteristic time 



of the matter advection due to the shear viscosity is t Vl 



j 



The time of the magnetic 



2 h B ^ 

field diffusion is tdiff = ^jif ? r l = 4^ = 6chv s . In the stationary state, the large-scale magnetic 
field in the accretion disk is determined by the equality t V i S = tdiff, what determines the ratio 



Br 

B 7 



avs_ 
av K 



a h 
a r 



<1 • 



(2.3) 



Here, v K = r£l K and j = rv K for a Keplerian disk. In a fully turbulent disk a matter is penetrating 
through magnetic field lines, almost without a field amplification. Note, that the field induced by 
the azimuthal disk currents has B z d ~ B r d [JTJl . 
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Figure 2: Sketch of the large-scale poloidal magnetic field threading a rotating turbulent accretion disk with 
a radiative outer boundary layer. The toroidal current flows mainly in the highly conductive radiative layers. 
The large-scale (average) field in the turbulent region is almost vertical. 



3. Turbulent disk with radiative outer zones 



Near the surface of the disk, in the region of low optical depth, the turbulent motion is sup- 
pressed by the radiative flux, similar to the suppression of the convection over the photospheres 
of stars with outer convective zones. The presence of the outer radiative layer does not affect the 
estimate of the characteristic time t v i sc of the matter advection in the accretion disk because it is 
determined by the main turbulent part of the disk. The time of the field diffusion, on the contrary, is 
significantly changed, because the electrical current is concentrated in the radiative highly conduc- 
tive regions, which generate the main part of the magnetic field. The structure of the magnetic field 
with outer radiative layers is shown schematically in Fig.2. Inside the turbulent disk the electrical 
current is negligibly small so that the magnetic field there is almost fully vertical, with B r <C B z , 
according to (^3|). In the outer radiative layer, the field diffusion is very small, so that matter ad- 
vection is leading to strong magnetic field amplification. We suppose, that in the stationary state 
the magnetic forces could support the optically thin regions against gravity. When the magnetic 
force balances the gravitational force in the outer optically thin part of the disk of surface density 
Z p h one finds the following relation takes place ^ 

r 2 ~ 2c ~ An' K } 

The surface density over the photosphere corresponds to a layer with effective optical depth close 
to 2/3 (e.g. [|l]]). We estimate the lower limit of the magnetic field strength, taking K es (instead of 
the effective opacity K e ff = ^/K es K a ). Writing K es ll p h = 2/3 , we obtain Z p h = 5/3 (g/cm 2 ) for the 
opacity of the Thomson scattering, K es = 0.4 cm 2 /g. The absorption opacity K a is much less than 



K es in the inner regions of a luminous accretion disk. Thus using in equation (pTl|), the above L p h, 



we estimate the lower bound on the large-scale magnetic field of a Keplerian accretion disk as 

c 2 1 _o 1 r M 

B z = W^^^ ^10 8 G— , %=- m= — . (3.2) 

The maximum magnetic field is reached when the outward magnetic force balances the gravi- 
tational force on the disk of surface mass density Z p h. In equilibrium, B z ~ y^Z^h. We find 
that B z in a Keplerian accretion disk is about 20 times less than its maximum possible value, for 
x = 10. a = 0.1, and m = 10. 
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4. Self-consistent numerical model 

Self-consistent models of the rotating accretion disks with a large-scale magnetic field requires 
solution the equations of magnetohydrodynamics. In presence of the radiative layer the strength of 
the magnetic field is large, and it may greatly exceed the strength of the seed field. The solution 
with a small field will not be stationary, and a transition to the strong field solution will take place. 
Therefore the strong field solution is the only stable stationary solution for a rotating accretion disk. 
The vertical structure of the disk with a large scale poloidal magnetic field was calculated in [|ll|], 
taking into account the turbulent viscosity and diffusivity, and the fact that the turbulence vanishes 
at the surface of the disk. The full system of equations was reduced to one vertical equation for the 
non-dimensional radial velocity u r , in the form 

4«2 d 2 ( d / d fl d / du r 

a P t79 P8^r PS 



+ « 2 /3 2 ^ (pg{u r - guo)) + ^ (^) + 3P&± = . (4.1) 

Here £ = z/h is a dimensionless height, u r = — v r / (occ s o), uo is a non-dimensional radial velocity 
in the non-magnetized disk [|l3|]. Coefficients of the turbulent viscosity v, and magnetic diffusivity 
T] are connected by the magnetic Prandtl number ^ ~1, v = ^T] = a ^ g(z) , where a is a 
constant, determining the turbulent viscosity jl3|]; j8 = c^ /v^ , where v^o = ^o/(47rpo) 1 ^ 2 is the 
midplane Alfven velocity, p = ^. The function g(z) accounts for the absence of turbulence in the 
surface layer of the disk [f|]. In the body of the disk g = 1, whereas at the surface of the disk, at 
say zs, g tends over a short distance to a very small value, effectively zero. The smooth function 

with a similar behavior is taken [[□]] in the form g(£)=^l — , with 8 <C 1. In the stationary 
state the boundary condition on the disk surface is u r = 0, and only one free parameter - magnetic 
Prandtl number remains in the problem. In a stationary disk vertical magnetic field has a unique 
value. The example of the radial velocity distribution for =1 is shown in Fig. 3. 



5. Jet collimation 

Magnetic collimation is connected with torsional oscillations of a cylinder with elongated 
magnetic field, see Fig.4. The stabilizing azimuthal magnetic field is created by torsional oscilla- 
tions. Approximate simplified model is developed [Q]. Ordinary differential equation is derived, 
and solved numerically, what gives a possibility to estimate quantitatively the range of parameters 
where jets may be stabilized by torsional oscillations. 

In non-dimensional variables T = (Ot. y = ■§- z = ar> = /?o = , differential 

equations have a form 

^=z,^ = i(l-Dsin 2 T),y(0) = l,z = 0atT = 0. (5.1) 
ax ax y 
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Figure 3: Distribution of the radial velocity over the thickness in the stationary accretion disk with a large 
scale poloidal magnetic field 




Figure 4: Jet confinement by magneto-torsional oscillations (qualitative picture) 

The system ( |57T1 ) has two non-dimensional parameters: D = 2 nKC, (^^) ' an( ^ an( ^ ^ e sec " 
ond one is taken equal to unity . Solution changes qualitatively with changing of the parameter D. 
At D < 2 there is no confinement, radius grows to infinity after several low-amplitude oscillations. 
At D = 2.1 radius is not growing to infinity, but is oscillating around some average value, forming 
complicated curves (Fig. 5).At D > 2.28 the radius goes to zero. On the edge of the cylinder the 
rotational velocity cannot exceed the light velocity. The analysis have shown [^] that for the sound 
velocity not exceeding c/2, the jet should contain baryons, density po exceeding about 30% of the 
total density of the jet. Development of chaos in this system, by construction of Poincare sections 
was investigated in [§]. 

6. Conclusions 

1. Disk field is amplified during disk accretion due to high conductivity in outer radiative 
layers. Stationary solution corresponds to j8 = 240 for Pr=l. 

2. Jets from accretion disk are magnetically collimated by large scale poloidal magnetic field 
by torsion oscillations, which may be regular or chaotic. 
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Figure 5: Time dependence of non-dimensional radius y (upper curve), and non-dimensional velocity z 
(lower curve), for D = 2.11 (left); for D = 2.25 (right), from [|] . 
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